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ABSTRACT 

The mean in-line force and the time dependent transverse 
force acting on a circular cylinder undergoing periodic trans- 
verse oscillations in an othem^ise steady flov; was measured 
for various amplitudes and frequencies of oscillation at several 
steady flow velocities. The experiments were carried out in a 
recirculating water tunnel operating as a closed channel at 
the test section. The mean in-line force has been expressed 
in terms of a mean drag coefficient and plotted as a func- 

tion of A/D and D/VT. The time dependent transverse force has 
been expressed in terms of the Fourier-averaged drag and inertia 
coefficients and ‘^ml and plotted as functions of the rela- 

tive amplitude A/D and the period parameter VT/D. The results 
have shown that the mean in-line force is significantly larger 
than that corresponding to steady flow about a non-oscillating 
cylinder and that energy may be transferred to the oscillating 
cylinder from the fluid at or near the vortex-shedding fre- 
quencies for the A/D values tested. 
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NOMENCLATURE 



A amplitude of cylinder oscillation 

Cdi mean in-line drag coefficient 

Fourier-averaged drag coefficient (Lift direction) 
Fourier-averaged inertia coefficient (Lift direction) 
D diameter of test cylinder 

F instantaneous total force acting on test cylinder 

T period of oscillation 

t time 

U instantaneous velocity 

maximum velocity in a cycle (U^=2 ttA/T) 

V mean velocity of the uniform flov; 

A/D dimensionless amplitude ratio 
D/VT dimensionless frequency parameter 
VT/D dimensionless period parameter 

Re Reynolds number (Re=VD/v) 

0 phase angle 

V fluid kinematic viscosity 

p fluid density 
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INTRODUCTION 



I . 

Elastic structures of one or more degrees of freedom can 
extract energy from the flow about them and can develop cata- 
strophic flow-induced vibrations. The understanding of this 
energy-extraction process is of paramount importance if one 
is either to eliminate or minimize it or to design the elastic 
structure such that it can withstand the oscillations under 
the contemplated environmental conditions. 

Even a superficial familiarity with the parameters involved 
in this type of phenomenon shows that the investigation of the 
flow- induced vibrations or the flow about oscillating bodies 
is more than an extension of the past studies on flow about 
bluff bodies at rest. In fact, the complexity of the problem, 
is increased by an order of magnitude. In view of this fact, 
it is no wonder that the past decade has produced either num- 
erous experimental data with the goal of obtaining ad-hoc 
solutions for specific problems, while directing little if 
any attention toward elucidation of the underlying mechanics 
and toward generalization of the results for future applica- 
tions; or numerous analytical models m.ost or all of which had 
nothing to do with the motion of the very medium which pro- 
vided the necessary energy to set the body in motion. In fact, 
the fluid mechanics of the flow-induced oscillations became 
so incidental to the phenomenon that the models dealt essenti- 
ally with black-box-induced oscillations. To be sure, these 
models heavily relied on experimental data partly to justify 
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their existence and partly to assign numerical values to 
numerous variable constants imbedded in them. The ability 
of these models to scale, i.e., to permit extrapolation, is 
uncertain. Thus, in the final analysis one is not quite 
sure whether one should use the experimental data within the 
range of their application or the curves fitted to them by 
the empirical models. Suffice it to say that the phenomenon 
is far from understood even for the idealized conditions 
encountered in the laboratory without the alarming consequences 
of the real ocean environment where fauna and flora, ever- 
changing ocean currents, and temperature gradients add further 
complications to an already complex problem. Surely, over 
design is not the answer but it may, in the next decade or so, 
be the best available tradeoff with failure. 

Mathematical models of flow-induced vibrations of bluff 
bodies and the response of circular cylinders to vortex shedding 
have been aptly described by Parkinson [1] and Currie, et al . 

[2] and will not be repeated here. 

It appears that among the various models considered so far, 
the "wake-oscillator" model of Hartlen and Currie [3] attracted 
more attention among the students of vibration analysis. We 
will not elaborate here on the attem.pts made by others to 
introduce one or more additional terms into the model origi- 
nally proposed by Hartlen and Currie for such naive attempts 
produced only more papers and unrealistically defined coeffici- 
ents than sound information toward the understanding of the 
flow-induced oscillations both in-line with and transverse to 
the ambient flow. 
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It suffices to conclude that despite the extensive atten- 
tion given to the study of the transverse oscillations of 
bluff bodies in air and water, there has been very little or 
no attempt to measure the forces acting on the bodies and to 
separate the time-dependent force into suitable components. 

It has, however, been firmly established that there is a band 
of frequencies in the vicinity of the Strouhal frequency in 
which the vortex shedding is synchronized or locked-in and 
that both the intensity of the vortex shedding and the span- 
wise correlation are enhanced. The synchronization or lock-in 
is manifested by a rapid decrease in inertial force and a 
rapid increase in the absolute value of the drag force. In 
other v;ords, the total force which is nearly in-phase with the 
motion before synchronization becomes nearly out of phase at 
and after synchronization. 

The purpose of the present investigation is the determina- 
tion of these in-phase and out-of-phase components of the time 
dependent force acting on cylinders undergoing transverse 
oscillations in a uniform stream and the correlation of the 
normalized force components with appropriate governing param- 
eters through the use of an equation proposed herein. Evi- 
dently, it is through the understanding of the forces acting 
on oscillating bodies that one would eventually be able to 
devise models to predict the strumjning or the hydroelastic 
behavior of rigid or elastic bodies immersed in fluid streams. 
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II. EXPERIMENTAL APPARATUS AND PROCEDURE 



A. EQUIPMENT 

1 . NPS Water Tunnel 

The experiments were performed in a recirculating v;ater 
tunnel (Figure 1) which had a capacity of approximately 500 
gallons. The galvanized test section, four inches wide, eight 
inches high and sixteen inches long, was closed on top with a 
removable plexiglass plate to eliminate the free surface 
effects. A small space adjacent to the side walls was provided 
to allow the passage of thin leaf arms connecting the cylinder, 
in the test section, to the driving hardware above the test 
section. A low RPM, high capacity, fourteen-inch-diameter- 
discharge centrifugal pump was used to circulate the fluid 
through the test section. The velocity of the fluid was regu- 
lated by a butterfly valve arrangement situated downstream of 
the test section. Velocities of 0.8 to 1.4 feet-per-second 
were obtained by adjusting the vanes of the butterfly valve. 

2 . Motor, Fly\^7hee'^ and Pivot Assembly 

The periodic motion in this experiment was obtained 
by the use of a small, variable speed, electric motor, a fly- 
wheel and pivot assembly (Figure 2) . The electric motor, 
mounted so as to isolate its vibrations from the remainder of 
the test apparatus, transmitted rotary motion through a flex- 
ible belt to a flywheel. A 56 inch long, rigid bar, pivoted 
at its center, transmitted the vertical component of motion 
of a bearing, attached to the flywheel, to a shaft connected 
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Figure 1. NPS water tunnel. 
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Fig. 2 Motor, Flywheel and Pivot Assembly 
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to the yoke assembly and cylinder. Horizontal motion was 
expended at both ends of the pivot rod by the free motion of 
bearings along the pivot arm. The amplitude, "A", of the 
vertical periodic motion was set by adjusting the radial posi- 
tion of the bearing attached to the flywheel. For these 
experiments, the amplitude ranged from 0.25 to 0.75 inches. 

The frequency of oscillation was regulated up to four cycles 
per second by adjustment of the variable speed motor. The 
system was designed and constructed to produce essentially 
sinusoidal oscillations, free of secondary oscillations. 

3 . Test Specimen and Yoke Assembly 

A one-inch diameter circular cylinder, constructed of 
aluminum tubing, was used as the test specimen. The cylinder 
was held in the test section by a yoke assembly which was 
connected to the pivot arm by means of a vertically constrained 
rod. The action of the pivot arm caused the cylinder to oscil- 
late vertically in the test section, transverse to the stream 
flow. The aluminum yoke assembly/force transducer (Figure 3) 
was instrumented with the strain gages and accelerometers 
required to monitor the forces and the acceleration felt by 
the cylinder. 

4 . Sensors and Sensor Support 

The mean velocity of the fluid was derived from the 
pressure differential of a pitot tube installed in the test 
section upstream of the test cylinder. One channel of a two 
channel recorder was used to continuously monitor the stream 
velocity during experiments. The calibration of the pitot 
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Fig, 3 Aluminum Yoke— Force Transducer 
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tube and transducer was accomplished by means of a water bottle 
and precisely measured plates which were used to produce a 
known water height differential. 

The instantaneous acceleration of the pivot arm was 
sensed by an accelerometer mounted on a yoke assembly. The 
acceleration signal was processed through a filter, to remove 
high frequency vibrations, and was recorded on one channel of 
a two channel recorder. 

Instantaneous forces were sensed by four piezo- 
resistive strain gages. A cantilever beam type arrangement, 
forming the top of the yoke assembly, as seen in figure 3, 
provided the mounting location for the strain gages. The 
forces sensed by this force transducer system were recorded 
on the two channel recorder. Calibration of the force trans- 
ducer system was accomplished by hanging loads at the mid- 
section of the cylinder, supported at the correct orientation 
so as to calibrate for vertical forces for the transverse or 
lift force measurements, and then for the horizontal, or in- 
line force measurements. 

B. PROCEDURE 

System checkout and sensor calibration preceded the 
experimental runs. Several test runs were made to verify the 
proper operation of the test system, sensors, processing and 
recording equipment, and to establish the test run procedures 
to be used. It was verified that the mechanical system pro- 
duced the required periodicity of oscillation for frequencies 
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belov/ about three cycles per second and amplitudes below one 
inch. It was decided that the system would not be driven at 
higher frequencies or larger amplitudes during the test runs 
because of the presence of excessive vibrational effects be- 
yond these limits. It was verified that the sensors and pro- 
cessing equipment had adequate sensitivity and phase reproduction 
to allow measurement of the required physical parameters. 

The first set of experiments has as an objective, the 
measurement of the mean fluid- induced force on the cylinder, in 
the direction of the stream flow, (in-line force) for various 
frequencies and amplitudes of cylinder oscillation transverse 
to the stream flow. It was verified by experiment that there 
was no in-line component of inertial force. 

The physical recording of the in-line force was accom- 
plished with the Honeywell recorder. Pitot differential pres- 
sure was recorded simultaneously to verify the magnitude and 
constancy of fluid velocity during the runs. Accelerometer 
output was recorded on the two channel recorder to provide 
data on oscillation period and phase angle. The flov; velocity 
was set at 0.84 feet per second, the amplitude at 0-25 inches, 
and the data runs were taken as the oscillation frequency was 
varied from zero to about three cycles per second. The charts 
were annotated with run data for later correlation. The pro- 
cedure was repeated for amplitudes of 0.5 and 0.75 inches and 
the entire series was repeated again at a flow velocity of 1.3 
feet per second. 



18 




h« 

4 

X 

r . 




♦ 

«7 JPUP ' 

" ■ .;;-• tap 

■**n4s^ 

»^ait 

1 

• » (v* ■• 

If Ud 

*.Mf> 

K 1 

» 4 « " 

I 

i 

Oltx 

1 

i A»* 



I 



Subsequently, the mean in-line force, the flow velocity, 
the oscillation period and amplitude, and the cylinder diam- 
eter were recorded from the charts. A mean drag coefficient 
C,. was manually calculated from this data for correlation 
with the frequency parameter D/VT and the nondimensionalized 
amplitude A/D. 

The second set of experiments had as an objective the 
measurement of the time dependent forces in the direction of 
cylinder m.otion, i.e., transverse to the direction of stream 
flow. Since, in this case, there was an inertial force in the 
direction of interest, it was necessary to adjust the experi- 
mental procedure in order to identify the instantaneous value 
of the inertial force. Separation of the inertial force from 
the total force measured with the cylinder submerged, v;as 
required to ascertain the resistive force produced by the 
fluid alone. 

The physical recording of the transverse force was accom- 
plished using one channel of the two channel recorder. Accel- 
eration was measured simultaneously on the other channel. The 
water level in the tunnel was alternately raised and lowered 
for each run in order to allow recording of the total "Wet" 
force and the inertial "Dry" force at a common frequency and 
amplitude of oscillation. The acceleration trace provided 
the reference for the time correlation of the wet and dry 
force data. Because of the necessity to properly correlate 
wet and dry phase information, extreme care was taken in veri- 
fying the zeros for the acceleration and force recordings. 
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since the time dependent force was of interest, the mean 
force and acceleration was used as the respective trace zeros . 
This value was found by electrically averaging the force and 
acceleration signals, and setting the resultant mean value to 
a convenient position on the recording trace. The flow vel- 
ocity was initially set at 0.84 feet per second, the amplitude 
at 0.25 inches, and runs were taken for various oscillation 
frequencies up to about three cycles per second. The procedure 
was repeated for amplitudes of 0.5 and 0.75 inches and then the 
entire series was repeated at a new flow velocity of 1.33 feet 
per second. Figure 4 shows the force and acceleration traces 
for a typical experimental run. 

From traces similar to figure 4, the transverse force was 
read and recorded from both the wet and dry run curves. Start- 
ing from the zero acceleration time (illustrated on figure 4) 
the wet and dry forces were recorded each 0.01 second, for at 
least three complete cycles. This data was then punched on IBM 
computer cards for evaluation. 

The computer program for calculating the transverse drag 
and inertia coefficients, given in Appendix A, was also used, 
along with visual inspection, to determine the goodness of the 
data and the accuracy of the data retrieval procedures. Com- 
puted values of the mean wet and dry forces, if significantly 
different from zero, indicated possible errors in the establish- 
ment of system zeros or the presence of secondary undesirable 
effects from the mechanical system. After completion of this 
data assessm.ent, the drag and inertia coefficients were 
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Fig, 4 Force and Acceleration Traces 
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calculated and recorded as a function of the period parameter 
VT/D for the three values of the nondimensionalized amplitude 
A/D investigated. 
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III. DATA ANALYSIS 



The numerical calculations as well as measurements in 
time-dependent flow yield the resultant force as a function of 
time for a given set of numerical values of the independent 
parameters. Thus it is not possible without a suitable hypo- 
thesis, to express the force both as a function of time and 
remaining independent parameters as one ordinarily would in a 
closed form solution. Such a working hypothesis is particu- 
larly necessary for the cable strumming problem since the 
results are to be incorporated into the dynamics of the cylin- 
der motion in the form of a forcing function. It should be 
stated at the outset that there is, at present, no generally 
accepted hypothesis to decompose the time dependent force into 
suitable components. As it will be seen shortly, even the 
existing hypotheses, such as the so-called Morison's equations, 
are not applicable to periodic flows with a non- zero mean 
velocity . 

Stokes, in a remarkable paper on the motion of pendulums, 
showed that the expression for the force on a sphere oscilla- 
tion in an unlimited viscous fluid consists of two terms, one 
involving the acceleration of the sphere and the other the 
velocity. This analysis shows that the inertia coefficient 
is modified because of viscosity and is augmented over the 
theoretical value valid for irrotational flow. The drag co- 
efficient associated with the velocity is modified because of 
acceleration, and its value is greater than it would be if the 
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sphere were moving with constant velocity. In general, the 
force experienced by a bluff body at a given time depends on 
the entire history of its acceleration as well as the instan- 
taneous velocity and acceleration. Thus, the drag coefficient 
in unsteady flow is not equal to that at the same instantaneous 
velocity in steady flow. Neither is the inertia coefficient 
equal to that found for unseparated potential flow. As yet a 
theoretical analysis of the problem for separated flow is 
difficult and much of the desired information must be obtained 
both experimentally and numerically. In this respect, the 
experimental studies of Morison and his co-workers [4] on the 
forces on piles due to the action of progressive waves have 
shed considerable light on the problem. The forces are divided 
into two parts, one due to the drag, as in the case of flow 
of constant velocity, and the other due to the acceleration 
or deceleration of the fluid. The concept necessitates the 
introduction of a drag coefficient and an inertia coeffici- 
ent in the expression for force. In particular if F is the 
force per unit length experienced by a cylinder, then 



F = 0 . 5 C^pD 



U 



D 



U + C pir . 
m 4 



dU 

dt 



( 1 ) 



where U and dU/dt represent respectively the undisturbed 
velocity and the acceleration of the fluid. 

On the basis of irrotational flow around the cylinder, 
should be equal to 2 (cylinder at rest, the fluid acceler- 
ating; otherwise = 1) , and one may suppose that the value 
of should be identical v/ith that applicable to a constant 
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velocity. However, numerous experiments show that this is 

not the case and that C. and C show considerable variations 

d m 

from those just cited above. Even though no one has suggested 
a better alternative, the use of the Morison ' s equation gave 
rise to a great deal of discussion on what values of the two 
coefficients should be used. Furthermore, the importance of 
the viscosity effect has remained in doubt since the experi- 
mental evidence published over the said period has been quite 
inconclusive . 

The drag and inertia coefficients obtained from a large 
number of field tests, as compiled by Wiegel [5], shov; exten- 
sive scatter whether they are plotted as a function of the 
Reynolds number or the so-called period parameter U^T/D. The 

reasons for the observed scatter of the coefficients C and 

m 

remained largely unknown. The scatter was attributed to 
several reasons or combinations thereof such as the irregu- 
larity of the ocean v-zaves, free surface effects, inadequacy 
of the average resistance coefficients to represent the actual 
variation of the nonlinear force, omission of some other 
important parameter which has not been incorporated into the 
analysis, the effect of ocean currents on separation, vortex 
formation, and hence on the forces acting on the cylinders, 
etc . 

The most systematic evaluation of the Fourier-averaged 
drag and inertia coefficients has been made by Keulegan and 
Carpenter [6] through measurements on submerged horizontal 
cylinders and plates in the node of a standing wave, applying 
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theoretically derived values for velocities and accelerations. 
Additional measurements have been made by Sarpkaya L7] of the 
in-line as well as transverse forces acting on cylinders and 
spheres in a sinusoidally oscillating fluid and it was found 
that the drag coefficient as well as the inertia coefficient 
for a strictly sinusoidally oscillating fluid (no mean vel- 
ocity) is a function of U^T/D and that the effect of the 
Reynolds number is rather secondary and certainly obscured 
by the excellent correlation of the data with the period param- 
eter U T/D. 
m 

On the basis of the above discussion, one would assume 
that Morison's equation would apply equally well to periodic 
flow with a mean velocity where u = V - U cos 6 and that C., 
and will have constant, time-invariant, Fourier or least- 

squares averages. This, in turn, implies that and 

are independent of the associated flow phenomena. There is, 
however, no a priori assurance in the principles of fluid 
mechanics of theory of models that this is, in fact, the case. 
Thus the effect of the combination of a uniform current and 
harmonic oscillations on the time-average and oscillatory 
forces acting on circular cylinders will have to be re-examined 
and the limits of application of the Morison's equation be 
delineated. 

It is a priori evident that both u = -U^cos6 and 
u = V - U^cos 9 yield the same acceleration du/dt. Thus, the 
force in-phase with the acceleration in Morison's equation 
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remains unaffected by the presence of the mean flow. The 
results presented herein show that this is not the case. 
Furthermore, the use of the Morison's equation as 




( 2 ) 



requires that the time-averaged drag force be calculated by 
increasing the force calculated from the steady flow by a 



that such an analysis appreciably underestimates the measured 
mean forces. It suffices to state that the fluid flow phen- 
omena for bluff bodies are significantly affected by the 
combination of currents and harmonic oscillations and that 
the results for steady currents alone and oscillations alone 
cannot be combined to yield reliable estimates of forces due 
to both acting together. 

The time-dependent forces per unit length in the present 
study are analyzed according to the following three-coeffici- 
ent equation: 



_ 2 

factor II + 0.5(U^/V) ]. The results presented herein show 



F = 0.5C,pDV^ + C pTT§- ^ (-U cosl^t) 
d m 4 dt m T 




(3) 



which may be written as 



0. 5PDV2 



F 



C , + C TT^(U T/D) (D/VT)^sin^t 
d ram'' T 



Cd (U^T/D) ^ (D/VT) ^ I cos|--t j cos^^t 



(4) 
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in which and are given by their Fourier averages as 




( 5 ) 



0 



and 



= -(3/4) y' (Fcos 0)d0/(U^^D) 
0 



( 6 ) 



Evidently, C^, C^, and are functions of VT/D and 
U^T/D or A/D. They may depend also on the Reynolds number 
which does not explicitly appear in the above expression be- 
cause of the assumptions made in the formulation of the basic 
force equation. 

In the foregoing, neither the coefficient is assumed 
to be equal to the steady-state drag coefficient for a uniform 
flow’ at the constant velocity V, nor C and C, are assumed to 



be identical to those obtained for a strictly harmonic oscil- 
lation. In fact, the results show that (steady) only 

for U =0, and and C are equal to those obtained for the 
m d m ^ 

harmonic oscillation only for VT/D=0. 

The equation proposed above is general enough to be appli- 
cable to both in-line and transverse oscillations. In the case 
of transverse oscillations, however, the mean net force in the 
direction of oscillation is zero, i.e., thus, one has 



m 



F (transverse force ) 
0. SpDV^L 



CmiTT^ (UmT/D (D/VT) ^sin^^t 



Cdi (UmT/D) ^ (D/VT) ^ |cos^t cos|^t 



(7) 
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In the foregoing, the inertia and drag coefficients have been 



denoted as C , and in order to distinguish them from those 

ml dl 

corresponding to in-line oscillations. The subscript "1" 
carries the meaning of "Lift" or force in the direction trans- 
verse to the stream. 

Ordinarily, for a perfectly sinusoidal oscillation of the 
cylinder, the coefficients and would be given by 

equations (5) and (6). However, when the oscillations are 
not perfectly harmonic, it is relatively more accurate to use 
the velocities and the accelerations encountered in the exper- 
ments rather than assuming them to be simple harmonic motions. 
It is with this objective in mind that the equations (5) and 
(6) were rewritten as 



as they would be obtained from equation (3) in the usual 
application of the Fourier analysis. Evidently, had the oscil- 
lations been perfectly harmonic the integrals appearing in the 
denominators of the above equations would have reduced to 



C 



2T 




( 8 ) 



ml 




and 




0 



(9) 



2tt 




( 10 ) 
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( 3 ) 




and 



2tt 

Jcos^ejcos o|d0 = I (11) 

0 



The equations (8) and (9) together v;ith equations (10) and 
(11) would have reduced to equations (5) and (6) . 

Since in the present study the periodic oscillations were 
not perfectly sinusoidal partly by design in order to obtain 
greater generality and flexibility in the experimentation and 
partly due to the constraints imposed in the design of the 
oscillation mechanism, it became necessary to incorporate into 
equations (8) and (9) the exact form of the oscillations 
encountered in the experiments. For this purpose the dry force, 
which is an exact representation of the oscillations , v;as norm- 
alized as 



f(6) = 



jF-dry (maximum)j 
and then the equations (8) and (9) v;ere revjritten as 



2ir 



'ml 



2T n / F-f(6)de 



piT 



and 



2d^U 

m 

f(0)d0 

0 

2it 

^2 q/f-I (0+7T/2)d0 



'dl 



2 2 

2fDTT A 



J^f (6 + ti/ 2) I f (6 + tt/ 2)| d0 



( 12 ) 



(13) 



(14) 
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It is easy to show that equations (13) and (14) reduce to 
equations (8) and (9) or to equations (5) or (6) for purely 
harmonic oscillations for which f(0) = sin 6. The advantage 
of the use of the equations (13) and (14) is rather obvious 
for all types of periodic oscillations. Furthermore, the 
independent evaluation of the denominators of equations (13) 
and (14) and their comparison with tt and 8/3 respectively (as 
they would have been equal to had the oscillations been har- 
monic) gave an indication of the deviation of the observed 
periodic oscillations from a purely harmonic oscillation. 

The force acting on the cylinder in the in-line direction 
due to the oscillations in the transverse direction is expressed 
in terms of a mean drag coefficient, denoted by given as 



^di 



Force in the in-line direction 



0 . 5 pLDV 



(15) 



Evidently, and depend on D/VT, U^T/D or 2 ttA/D, 

and possibly on the Reynolds niimber. 

The experimental data are analyzed using the computer pro- 
gram given in Appendix A written according to the equations 
(13) , (14) , and (15) and are plotted in terms of A/D and 

either D/VT or VT/D. 
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IV. DISCUSSION OF RESULTS 



The results v/ill be discussed in two parts. The first 
will be the average in-line force acting on the cylinder 
undergoing forced periodic oscillations in the transverse 
direction. The second will be the time dependent transverse, 
force . 

Evidently, the average in-line force is coupled with 
secondary oscillations due to vortex shedding. Hov;ever, 
such oscillations are rather small in both steady and peri- 
odic flows and certainly not larger than about seven percent 
of tlie average force. It is for this reason that only the 
average of the in-line force acting on the oscillating cylin- 
der is presented here. 

Figures 5, 6, and 7 show the variation of the normalized 
in-line force as a function of D/VT for A/D =0.25, 

0.50 and 0.75 respectively. Each figure represents the data 
obtained with two velocities, namely, V=0.84 and V=1.3. In 
normalized form these velocities correspond, for the one 
inch cylinder used, to the Reynolds numbers Re=VD/v = 7000, 
and Re = 10,833. 

Evidently, the in-line force increases with A/D since 
the cylinder, undergoing transverse oscillations, presents 
a larger apparent-pro jected area to the mean flow. This, 
however, is only part of the explanation. In addition, the 
vortex grov/th and motion are affected by the oscillation of 
the cylinder which in turn affect the in-line and ti'ansverse 
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forces acting on the cylinder. This is evidenced by the fact 
that the in-line force for a given A/D increases at first, 
reaches a maximum, and then decreases as D/VT increases. A 
simple minded calculation based on the steady flow drag co- 
efficient for a stationary cylinder and the apparent projected 
area for the in-line force F, which may be written as 



'di 



JspV DL 



+ 2A/D) 



yields values which are almost equal to the maximum values 
given in figures 5, 6, and 7. It should be noted, hov/ever, 
that the phenomenon is far more complex, and that such a 
simple minded procedure should not generally be used, even 
though the results are surprisingly good. 

For the purposes of comparison, the figures 5, 6, and 7 
are combined in figure 8 by drawing mean lines through the 
data points. Figure 8 shov;s that the in-line force coeffici- 
ent reaches its maximum at D/VT between 0.18 and 0.20. Ordin- 
arily, the Strouhal number for a stationary cylinder would be 
0.22 for the Reynolds numbers cited previously, and one would 
expect that the forces acting on the cylinder will undergo 
dramatic changes as the vortex shedding frequency given by 
the Strouhal number coincides with the frequency of the cylin- 
der oscillations. The present results show that such a syn- 
chronization takes place at a frequency slightly lower than 
the Strouhal frequency. 

Also shown in figure 8 is the normalized amplitude of 
the oscillations in the in-line force for A/D=0.5. As noted 
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earlier the oscillations are quite negligible relative to the 
mean force and certainly under 7 per-cent. It should be noted 
that the amplitude of the oscillations, like the average force 
show an almost sudden increase in the vicinity of D/VT nearly 
equal 0.19 and remain at that value for larger values of D/VT. 
The occurrence of synchronization as well as the increase of 
the amplitude of oscillations are shown m.ost dramatically in 
figure 9. This figure was obtained by setting the free stream 
velocity at 0.84 feet per second and the A/D ratio equal to 
0.5. Then, beginning with the case of the non-oscillating 
cylinder, the frequency of the oscillations was gradually 
increased up to about four cycles per second and the resulting 
in-line force was continuously recorded. The figure shows 
that the in-line force increases rapidly but with very little 
oscillations superimposed on it. As soon as the state of 
synchronization is reached, the amplitude as well as the fre- 
quency of the force oscillations increases. 

From an engineering point of view the significance of the 
magnitude of the in-line force is that a cylinder or cable 
excited by the flow to oscillate in the transverse direction 
may be subjected to in-line forces several times larger than 
assumed in its design. Furthermore, the deflections caused 
by the in-line force on sufficiently flexible cylinders tend 
to couple with transverse ocsillations and not only affect 
the magnitude of the transverse oscillations but also the 
path of the cylinder motion. Thus it is not uncommon to see 
heat exchanger pipes or chimneys exhibit oscillation patterns 
in the plane normal to their axis. 
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The time dependent transverse force is described, as 

noted earlier, in terms of a drag coefficient and an 

inertia coefficient given by 

2 2 I 

Transverse Force=C , ptiD d (-U cos6) - C,,pDU cos 6 cos 0 

ml-j — m dl TT m J 



Figures 10, 11, and 12 show and as a function of 

VT/D for A/D - 0.25, 0.50 and 0.75 respectively. These co- 
efficients were obtained without the term in the general- 
ized Morison equation (see equation 3) since the mean of the 
transverse force is zero. 

It is seen from these figures that important variations 
in and occur particularly in the vicinity of the 

Strouhal frequency (here VT/D=4 . 5 to 7) where the natural 
eddy-shedding at the Strouhal frequency is both enhanced and 
correlated by the oscillations. 

The inertia coefficient or the normalized in-phase com- 
ponent of the transverse force undergoes a rapid drop as the 
frequency of the oscillation approaches the Strouhal frequency 
from both the upper and lower limits. In other words, syn- 
chronization or lock-in is manifested by a rapid decrease 
in inertial force and a rapid increase in the absolute value 
of the drag force. 

This fact, which has not been recognized before, shows 
that the lock-in phenomenon is a phase transformer. The 
total force which is nearly in phase with the motion before 
synchronization becomes nearly out of phase at or after syn- 
chronization. It should be noted in passing that the success 
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of many empirical models dealing with this type of oscillation 
comes partly from the manipulation of the phase angle near 
synchronization. It is now apparent that the fluid force to 
be used in the equations expressing the motion of a cylinder 
or cable should be given by the data presented herein. Such 
data take care, not only of the variation of the phase angle, 
but also the amplitude of the transverse force as a function 
of the normalized frequency. 

Figure 13 depicts an example of the occurrence of syn- 
chronization as the period of oscillation is gradually de- 
creased. The upper trace shows the phenomenon as frequency 
is increased from the non-oscillating case to just beyond 
the synchronization frequency. The lower trace provides an 
exploded view of the phenomenon near the synchronization fre- 
quency. The transverse force sharply decreases as soon as 
the point of synchronization is reached. Beyond that point, 
the changes in the transverse force with frequency are quite 
small. The phenomenon is reversible and one would obtain a 
figure similar to figure 13 if one gradually decreased the 
frequency. The possibility of the occurrence of an hysteresis 
has not been investigated. 

The drag coefficient or the normalized out-of-phase 

component of the total instantaneous transverse force given 
in figures 10, 11, and 12 show that becomes negative for 

VT/D values between approximately 4.5 and 7 (i.e. for norm- 
alized frequencies between 0.14 and 0.22). Outside this 
range the drag is positive, thus in the opposite direction 
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to the motion of the cylinder. Within the range of VT/D values 
cited above, the drag force is in phase with the direction of 
motion of the cylinder and helps to magnify the oscillations 
rather than damp them out. For this reason, the region in 
which is negative is sometimes referred to as the negative 

damping region. The fact of the matter is that this is not 
damping in the proper use of the v/ord but rather an energy 
transfer from the fluid to the cylinder via the mechanism of 
synchronization. The values of VT/D at which changes its 

sign depend on A/D as seen in figures 10, 11, and 12 even 
though the negative region roughly lies within the VT/D 

values of 4.5 and 7. The envelopes of the two zero crossings 
determine the region of self excited oscillations. A precise 
determination of such an envelope and its dependence on the 
Reynolds number will require experiments with additional A/D 
values and other Reynolds numbers. For purposes of the present 
investigation, it suffices to note that the maximum absolute 
value of in the synchronization range decreases rapidly 

as A/D increases. Field studies have shown that synchroniza- 
tion does not occur for relative amplitudes larger than unity. 
The trend of the present data is in conformity with such 
observations . 

Finally, an unexpected and previously unknown observation 

in connection with the variation of C,, will be described. 

dl 

For normalized frequencies (D/VT) smaller than about 0.14 
the data yield positive drag coefficients. Between 0.14 and 
0.22, the synchronization takes place and the drag coefficient 



37 









I 



« i >. 



■ i. 



t 

r 



♦ 



is negative as noted above. Ordinarily one would have expected 
that the drag coefficient will remain positive and continue 
to increase with increasing frequencies beyond D/VT = 0.22 and 
eventually reach a value v;hich would be identical to that 
obtained by oscillating the cylinder in a fluid otherwise at 
rest. However, an interesting phenomenon takes place at fre- 
quencies between approximately 0.22 and 0.27, depending on 
the A/D ratio. For example for A/D=0.75, increases 

rapidly at D/VT = 0.26 and then decreases sharply to nearly 
zero. At first it was suspected that this might be due to an 
experimental error. However, the repeatability of the experi- 
ments and the observation of the same phenomenon at other 
relative amplitudes and velocities have shown that there is 
indeed a dramatic change in at D/VT = 0.26 for A/D = 0.75, 

at D/VT =0.24 for A/D = 0.5 and at D/VT= 0.31 for A/D = 0.25. 
It is further noted from the data presented herein that, 
particularly for A/D = 0.25 and 0.5, becomes once more 

negative in a narrower range of D/VT values (D/VT from 0.345 
to 0.4 5 for A/D = 0.25) and shov;s the existence of a second 
region of synchronization. The D/VT value at which 
acquires its second minimum value in the case of A/D = 0.25 
is almost exactly twice that of the first minimum. The 
occurrence of this second region of synchronization at higher 
frequencies depends on the A/D value. The fact that is 
demonstrated here is that there is not a single region of 
synchronization and that there is at least one and possibly 
more regions of frequency in which synchronization can occur. 
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The narrowness of the regions of frequency in the second 
region of synchronization makes it rather difficult to observe 
the phenomenon. In fact one may easily miss such a region by 
simply not taking smaller increments in frequency. It suf- 
fices to say that a cylinder may be excited first at fre- 
quencies near the Strouhal frequency and then at the multiples 
of the Strouhal frequency. However, the largest energy trans- 
fer from the fluid to the cylinder occurs in the first syn- 
chronization region near the Strouhal frequency. 
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V. CONCLUSIONS 



The experimental investigation of the transverse oscilla- 
tions of a cylinder in a flow with an ambient mean velocity 

V has yielded the force coefficients C^., C , and C,, and 
^ di ml dl 

has shown that; 

a. The mean flow has significant effects on the force 
transfer coefficients and that the result of experiments with 
harmonic oscillations in a fluid otherwise at rest are not 
applicable to the transverse oscillations of a cylinder in a 
uniform flow; 

b. It is possible to excite transverse oscillations for 
A/D smaller than about unity. In a region of D/VT in the 
vicinity of the Strouhal frequency, the out-of-phase compo- 
nent of the total force becomes negative and some energy is 
transferred from the fluid to the cylinder. The rate at 
which energy is transferred decreases with increasing rela- 
tive amplitudes. Furthermore there is at least one, and 
possibly m.ore, narrower bands of frequencies at which syn- 
chronization occurs. 

c. Transverse oscillations give rise to an increa^ 
drag force in the direction of the mean flow. This force 
depends on A/D as well as on D/VT and reaches a maximum at 
about D/VT = 0.18. This value corresponds to a normalized 
frequency slightly belov; the Strouhal frequency and is well 
within the synchronization region. From a practical point of 
view this is a matter of concern for sound design of cables 



40 




\ 



;*• 





I 






and other structures which may be subjected to transverse 
oscillations? 

d. A meaningful dynamic analysis of the vortex excited 
transverse oscillations may be carried out only through the 
use of the force transfer coefficients presented herein. 
However, it is necessary that additional data be obtained at 
other amplitudes and Reynolds numbers. Furthermore, additional 
experiments may also have to consider the roughness of the 
oscillating structure. 
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Mean in-line drae coefficient versus D/VT for A/D 
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in-line drag coefficient versus D/VT for A/D = 0.75 
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Comparison curves: Mean in-line dreq coefficient versus D/VT 
for A/D = 0.25, 0.50, 0.75 
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Phenoncnological Dorrionstration : In-line force versus tine while 

increasing oscillation frequency, (A/D = 0.50, V = 0.84, D/VT = ^ 
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Fig, 10 Fourior-avcrdtjecl drag and inertia coefficients C ,, and. C versus VT/P for !^./D = 0.25. 
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Fig. n Fourier-averaged drag and inertia coefficients C.-, and C , versus VT/D for A/D 
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Fig. 12 Fourier-averaaed drag and inertia coefficients C,, and C •, versus VT/D for A/D = 0.75. 
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